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Abstract 

Through an extension of the ultradiscretization for the optimal velocity (OV) model, we introduce 
an ultradiscretizable traffic flow model, which is a hybrid of the OV and the slow-to-start (s2s) 
models. Its ultradiscrete limit gives a generalization of a special case of the ultradiscrete OV (uOV) 
model recently proposed by Takahashi and Matsukidaira. A phase transition from free to jam 
phases as well as the existence of multiple metastable states are observed in numerically obtained 
fundamental diagrams for cellular automata (CA), which are special cases of the ultradiscrete limit 
of the hybrid model. 
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I. INTRODUCTION 



Studies on microscopic models for vehicle traffic provided a good point of view on the 
phase transition from free to congested traffic flow. Related self-driven many-particle systems 
have attracted considerable interests not only from engineers but also from physicists 
Among such models, the optimal velocity model 3|], which successfully shows a formation 
of "phantom traffic jams" in the high-density regime, is a car-following model describing an 
adaptation to the optimal velocity that depends on the distance from the vehicle ahead. 

Whereas the OV model consists of ordinary differential equations (ODE), cellular au- 
tomata (CA) such as the Nagel-Schreckenberg model 1^, the elementary CA of Rule 184 
(ECA184) 5], the Fukui-Ishibashi (FI) model ^ and the slow-to-start (s2s) model [7j are 
extensively used in analyses of traffic flow. Recently, Takahashi and Matsukidaira proposed 
a discrete OV (dOV) model, which enables an ultradiscretization of the OV model [sl. The 
resultant ultradiscrete OV (uOV) model includes both the ECA184 and the FI model as its 
special cases. However, the s2s effect remains to be included in their ultradiscretization. The 
aim of this letter is to present an ultradiscretizable hybrid of the OV and the s2s models. 

II. THE OV MODEL AND THE S2S EFFECT 

Imagine many cars running in one direction on a single-lane highway. Let Xk{t) denote 
the position of the k-th car at time t. No overtaking is assumed so that Xk{t) < Xk+i{t) 
holds for arbitrary time t. The time-evolution of the OV model [s!] is given by 

dVk{t) 1 



dt 



-(vopt{AXk{t)) -Vk{t)), (1) 



where Vk ■= ^ and Axk := Xk+i — Xk are the velocity of the A;-th car and the interval 
between the cars k and /c + 1, respectively. A function fopt and a constant to represent an 
optimal velocity and sensitivity of drivers, or the delay of drivers' response, in other words. 

Since the current velocity and the current interval between the car ahead determine 
the acceleration through the time-evolution and the optimal velocity, we classify the OV 
model ([T]) as the acceleration-control type (aOV). On the other hand, the OV model of the 
velocity-control type (vOV) was proposed in earlier studies of the car-following models jo], 

Vk{t)=v,^,{Axk{t-to)). (2) 



Replacement of t in the above equation ([2]) with t + to and the Taylor series of Vk{t + to) 
yield 



Vn 



-,{Axk{t)) = Vk{t + to) = Vk{t) + — + + 



.optV -V" ■ "U, ^.v^, . dt ' 2 dt2 

which is rewritten as 



dt 2 dt2 



toH = ^(vopt {Ax k{t)) -Vk(t)). 

to V / 



Thus we note that the aOV model ([T]) is given by neglection of the higher derivatives in the 
Taylor series of the vOV model ([2]). Though the aOV model is more common in the studies 
on vehicle traffic, we shall concentrate on an ultradiscretizable hybrid of the vOV and the 
s2s models. Thus we call the vOV model ([2]) simply as the OV model, hereafter. 

Note that the input to the OV function Vopt{x) in the OV model ([2]) is the headway at 
a single point of time t — to that is prior to the present time t. Thus we may say that the 
OV model describes, in a sense, "reckless" drivers since the model pays no attention to the 
headway between the time t — to and the present time t. On the other hand, "cautious" 
drivers governed by the s2s model Tj keep watching and require enough length of headway 
to go on for a certain period of time before they restart their cars. The contrast between 
the two models suggests the idea that the s2s effect and the OV model can be brought 
together by appropriately choosing an effective distance AeffXfc(t) containing information on 
the headway for a certain period of time going back from the present as an input to the OV 
function Vopt{x). We shall see this idea works in what follows. 

What is crucial in the ultradiscretization of the aOV model [s^ is the choice of the OV 
function, 

where Vo,Xq and 6x are positive constants. In terms of the auxiliary functions, 

v.p.{x) := vo^^ (4) 
Jopt(a:) :=5xlog(l + e(^-^°)/^^) (5) 

the OV function ([3]) is expressed as 

Vopt{x) = Vopt{x) - Voptix = 0). 
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A naive discretization of the auxiliary function (jl]), 

5opt(a;) - £opt(a; - vo5t) 



6t 



introduces the OV function for the discrete OV (dOV) model, 



vtpti^) = vt^^ix) - v^^^{x = 0) 



^opt 



6x 



log 



]_ _j_ q{x-xo)/5x /]__!_ Q{x-xo-voSt)/5x 
1 _|_ Q-xo/Sx / ]^ _|_ Q-{xo+vo5t)/Sx 



(6) 



which is found to be ultradiscretizable. 



Let := Xk{t = n5t) and := {x^'^^ -xl)/6t where n{= 0, 1, 2 ■ • ■ ) and 6t{> 0) are the 
integral time and the discrete time-step, respectively. Employing the effective distance as 



A^ffX^ :=5xlog(^ 



n'=0 



no + 1 



where hq := to/ St, we extend the OV model ([2]) in a time-discretized form as 



(7) 



n „,d /Ad 



= ^opt 



(8) 



which is equivalent to 



X 



n+l n 



xl + 5x< lo£ 



!^ p-{Ax"-" -xo)/5x. _i 



E 

n'=0 



no + 1 



log(l +e-"o/^") 



-log 



"0 g-(Ax^-" -zo-t>o<5t)/<5x ^ -1 



n'=0 



no + 1 



+ log(l + e 



-(xo+?;o<5t)/(5a; 



It is straightforward to confirm that the continuum limit (5t — ^ of the above discrete s2s-0V 
(ds2s-0V) model (jHD reduces to the integral-differential equation which we call the s2s-0V 
model. 



dt 



Jopt 



{A,sXk{t)) =Vo(l + - r° e-(A-.(t-t')-o)/fe^^A ' (1 + ^x0/5x^-1^ 



where the corresponding effective distance is given by 



A,sXk := 5xlog(- / " e-^^k(t-t')/Sx^^f\ \ 
^^0 Jo ^ 



We shall see that the s2s effect is indeed built into the OV model in the ultradiscrete limit 
of the ds2s-0V model. 
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III. ULTRADISCRETIZATION 



Ultradiscretization lOj is a scheme for getting a piecewise-linear equation from a differ- 



ence equation via the hmit formula 

hm (5x(e^/''" + e^/^" + ■•■)= max(A, 5, • ■ • )• 

In order to go forward to the uhradiscretization of the ds2s-0V model (IHl), it will be a good 
choice for us to begin with the ultradiscrete limit 5x +0 of the auxiliary function ([5]): 



x"p^(x) := lim Xopt(a;) = max(0, x — Xq). (10) 



In the same way to make the OV function for the dOV model ([6]) from the auxiliary func- 
tion ([S]), we obtain the OV function for the uOV model as 

<pt(2;) = ^pt{^) - vlp^{x = 0) = max(^0, ^ - max(^0, ^ - wo), (H) 

where "yoptl^) (5^opt(^) ~ 5;"p^.(a; — fo5t))/5t. The effective distance ([7]), on the other hand, 
is ultradiscretized in the same manner: 

A^gx^ := lim A^ffX^ = -m;^x(-Axr"') = rnWAx—'). (12) 

ax— >+0 n'=0 n'=0 

Thus we obtain an ultradiscrete equation 



^" = <pt(A>^), (13) 



which is equivalent to 

1 + max^O, niin(Ax^~"') - Xoj - max^O, niin(Ax^""') - xq - voSt^ , 



^ri+l _ n 



as the ultradiscrete limit of the ds2s-0V model ([8]). We name it the ultradiscrete s2s-0V 
(us2s-0V) model. When the monitoring period no is fixed at zero, the us2s-0V model 
reduces to a special case of the uOV model Jsl. As we can see from eqs. f fTTj) . f|T2l) and 
(IT^ . the velocity is determined by the optimal velocity for the minimum headway in 
the period between n — hq and n. Thus cars will not restart nor accelerate, unless enough 
clearance goes on for a certain period of time. On the other hand, cars immediately stop 
or slow down when their headways become too small to keep their velocities. The s2s effect 
and a "cautious" manner of driving are built into the uOV model in this way. 



Now let us see how a CA comes out from the us2s-0V modeL Let xq be the discretization 
step of the headway Ax^, or equivalently, the size of the unit cell of the CA. Then with 
no loss of generality, we may set Xq = 1. Assume that the number of vacant cells between 
the cars k and k + 1, Ax"^ := Ax^ — xq, must be non-negative, Ax^ > 0, which prohibits 
car-crash. Then the us2s-0V model ffT3ll reduces to 



,1 ' "0 



^ +min(^mm(Ax^ "'),t;o(5tj. (14) 

Fixing vo5t at an integer, we call this model the s2s-0V cellular automaton (CA). The 
s2s-0V CA reduces to the FI model 6| when no = and to the ECA184 jsl when no = 
and v^dt = 1(= xq). The s2s model 7|] also comes out from the s2s-0V CA by choosing 
no = 1 and Vo6t = 1(= Xq). Thus the s2s-0V CA is regarded as a hybrid of the FI model 
and an extended s2s model. 



IV. NUMERICAL EXPERIMENTS 



We shall numerically investigate the s2s-0V CA ([HI). Throughout this section, the length 
of the circuit L is fixed at L = 100 and the periodic boundary condition is assumed as well 
so that x^ + L is identified with x^. 

Spatio-temporal patterns showing trajectories of each vehicle are given in fig. [TJ We 
choose the parameters and initial conditions so that jams appear in the trajectories. The 
two figures in the top share the same monitoring period no = 2 but their maximum velocities 
are different. The top left trajectories show that the velocities of the vehicles are zero or one, 
which is less than or equal to its maximum velocity vo6t = 1. In the top right trajectories 
whose maximum velocity VQ6t = 3, on the other hand, the velocities of the vehicles read 
zero, one, two and three. Thus we notice that the vehicles driven by the s2s-0V CA can run 
at any allowed integral velocity which is less than or equal to its maximum velocity v^St. 

The other two figures in the bottom in fig. [T] share the same maximum velocity voSt = 2, 
but their monitoring periods are different. As is observed in the bottom two figures, the 
longer the monitoring period is, the longer it takes for the cars to get out of the traffic 
jam. The jam front is observed to propagate against the stream of vehicles at constant 
velocity (^no+i)st ' since cars have to wait no + 1 time-steps to restart after their preceding 
cars restarted, as is depicted in fig. [2l 
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FIG. 1: The spatio-temporal patterns of the s2s-0V CA. For ah four patterns, the number of 
cars K is fixed at = 30 The maximum velocities vo6t and the monitoring periods no for these 
patterns are (top left) vo5t = 1, uq = 2, (top right) vo6t = 3, no = 2, (bottom left) vo5t = 2, no = 1 
and (bottom right) vo5t = 2, uq = 3, respectively. 



Fig. [3] shows fundamental diagrams giving the relation between the vehicle flow 



Q:-- 



{rii - no + l)L 



EE 



6t 



k=l n=no 

which is equivalent to the total momentum of vehicles per unit length, and the vehicle density 
p := j^, where K is the number of vehicles. The fundamental diagrams clearly show phase 
transitions from free to jam phases as well as metastable states, which are also observed in 
empirical flow-density relations . It is remarkable that the fundamental diagrams have 
multiple metastable branches. This feature is similar to that reported by Nishinari et al. 11 1 
By observation, we note that each fundamental diagram has v^dt metastable branches and a 
jamming line. The branches and the jamming line correspond to integral velocities that are 
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FIG. 2: Backward propagation of the jam front at constant velocity (n^j^x'^i^t 
vo5t = 2, no = 3 and xq = 1. 



J for the case 



less than or equal to the maximum velocity VoSt. Let us confirm it with fig. [3l The top two 
figures share the same monitoring period tiq = 3, but their maximum velocities are different. 
The top left diagram corresponding to VQ6t = 2 has three branches. This number equals 
to that of all the integral velocities, two, one and zero, as is depicted in the diagram. The 
number of the metastable branches in the top right diagram as well as those of the bottom 
two are explained in the same manner. This observation also suggests that the monitoring 
period is irrelevant to the number of metastable branches. 

All the end points of the branches as well as the jamming line are on the line p + Q{= 
p + Q^) = 1. This is because the density at the end point is the maximum density Pmax(^) 
that allows the velocity of the slowest car to be v6t. The maximum density Pmax(^^) is 
determined by 

vdt + Xo 

Since all the cars flow at the velocity v6t when p = pmax{v)i the corresponding flow is given 
by (5(Pmax) = ^^Pmax- Thus the relation p^ax + <5(Pmax)|J = 1 holds. 

The free line is a branch whose inclination equals to the maximal velocity v^dt. Any other 
metastable branch and the jamming line branch out from the free line. By observation, the 
density of the branch point of the branch corresponding to the velocity v5t reads 

_ 

[v^dt — v5t)no + v^dt + xq 
This observation is explained as follows. Suppose one car, say the car /c, runs at the velocity 
V and the other K — 1 cars run at the maximum velocity Vq. At the moment the /c-th car 
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Density Density 

FIG. 3: The fundamental diagrams of the s2s-0V CA. The flows Q arc computed by averaging 
over the time period 800 < n < 1000. The maximum velocities vqSI and the monitoring periods 
no for these patterns are (top left) v^dt = 2, no = 3, (top right) v^St = 4, no = 3, (bottom left) 
v^St = 3, no = 2 and (bottom right) v^St = 3, no = 4, respectively. The inclination of the free line 
equals to the maximum velocity VQ5t. The jamming line has a negative inclination. 

slows down to the headway between the cars k and /c + l is v5t-\-XQ. Since it takes at least 
no + 1 time-steps for the car k to speed up to vq, the headway between the cars k and k-\-l 
expands up to if = (fo^t — v5t){nQ + 1) + v5t + xq = xq/ p]^ > vo6t by the time the k-th car 
speeds up to Vq. If all the cars can obtain the headway H, slow cars running at the velocity 
V disappear in the end. Thus the density at the branch point of the branch corresponding to 
the velocity v5t is given by pb — xq/H. Note that the density at the branch point becomes 
smaller as the monitoring period becomes larger. 
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V. CONCLUDING REMARKS 



Through an extension of the ultradiscretization for the OV model [8j, we introduced the 
ds2s-0V ([HD and s2s-0V Qj models as ultradiscretizable traffic flow models. The model is 
a hybrid of the OV [3] and the s2s models whose ultradiscrete limit gives a generalization 
of a special case of the uOV model by Takahashi and Matsukidaira 8| . The phase transition 
from free to jam phases as well as the existence of multiple metastable states were observed 
in numerically obtained fundamental diagrams for the s2s-0V CA (1141) . which are special 
cases of the us2s-0V model f|T3|) . 

Detailed studies on the properties of the hybrid models ([H]), ©, (H^ and (HM such as 
exact solutions, comparison with other traffic flow models as well as empirical data remain 
to be investigated. 

Acknowledgments 

The authors are grateful to D. Takahashi, J. Matsukidaira, A. Tomoeda, D. Yanagisawa 
and R. Nishi for their valuable comments at the spring meeting of JSIAM in March, 2009. 



[1] D. Chowdhury, L. Santen and A. Schadschneider, Statistical physics of vehicular traffic and 

some related systems, Phys. Rep. 329 (2000), 199-329. 
[2] D. Helbing, Traffic and related self-driven many-particle systems. Rev. Mod. Phys., 73 (2001), 

1067-1141. 

[3] M. Bando, K. Hasebe, A. Nakayama, A. Shibata and Y. Sugiyama, Dynamical model of traffic 
congestion and numerical simulation, Phys. Rev. E, 51 (1995), 1035-1042. 

[4] K. Nagel and M. Schreckenberg, A cellular automaton model for freeway traffic, J. Physique 
I, 2 (1992), 2221-2229. 

[5] S. Wolfram, Theory and applications of cellular automata. World Scientific, Singapore, 1986. 
[6] M. Fukui and Y. Ishibashi, Traffic flow in ID cellular automaton model including cars moving 

with high speed, J. Phys. Soc. Jpn. 65 (1996), 1868-1870. 
[7] M. Takayasu and H. Takayasu, 1// noise in a traffic model. Fractals 1 (1993), 860-866. 



10 



[8] D. Takahashi and J. Matsukidaira, On a discrete optimal velocity model and its continuous 

and ultradiscrete relatives, JSIAM Letters 1 (2009), 1-4. 
[9] G. F. Newell, Nonlinear effects in the dynamics of car following, Oper. Res. 9 (1961), 209-229. 
[10] T. Tokihiro, D. Takahashi, J. Matsukidaira and J. Satsuma, From soliton equations to inte- 
grable cellular automata through a limiting procedure, Phys. Rev. Lett. 76 (1996), 3247-3250. 
[11] K. Nishinari, M. Fukui and A. Schadschneider, A stochastic cellular automaton model for 
traffic flow with multiple metastable states, J. Phys. A. Math. Gen. 37 (2004), 3101-3110. 



11 



